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Abstract. We present a detailed analysis of tlie medium redshift (z = 0.2906) galaxy cluster Z3146 which is part of the 
ongoing ARCRAIDER project, a systematic search for gravitational arcs in massive clusters of galaxies. The analysis of Z3 146 
is based on deep optical wide field observations in the B, V and R bands obtained with the WFI@ES02.2m, and shallow archival 
WFPC2@HST taken with the F606W filter, which are used for strong as well as weak lensing analyses. Additionally we have 
used publicly available XMM/Newton observations for a detailed X-ray analysis of Z3146. Both methods, lensing and X-ray, 
were used to determine the dynamical state and to estimate the total mass. We also identified four gravitational arc candidates. 
We find this cluster to be in a relaxed state, which is confirmed by a large cooling flow with nominal ~ 1 600 Mq per year, regular 
galaxy density and light distributions and a regular shape of the weak lensing mass reconstruction. The mass content derived 
with the different methods agrees well within 25% at = 1661 h^^^kpc indicating a velocity dispersion of cr,, = 869;^|j3 km/s. 
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1. Introduction 

Galaxy clusters are the largest bound structures in the uni- 
verse and therefore excellent cosmological probes. In partic- 
ular, large samples of clusters allow a statistical study of their 
physical properties. Such samples need clear selection criteria, 
e.g. selection by mass. Due to the tight relation between the X- 
ray luminosity Lx and the mass (Schindler 1999 Reiprich & 



|Bohringer 1999 1 X-ray surveys provide an excellent basis to 
select the most massive systems for lensing studies (IL uppino] 
[etaLj [T 9991 [Smith et al.[ [200T] |2005) [Bardeau et al. , 2005^ 
The combination of lensing and X-ray studies allows us to get 
important insights into galaxy clusters, as it offers the pos- 
sibility to obtain physical properties of the cluster members, 
the Intra Cluster Medium (ICM) and the determination of the 
cluster gravitational mass and its distribution with independent 
methods. 

Send offprint requests to: W. Kausch 

* Based on observations made with the NASA/ESA Hubble Space 
Telescope, obtained from the data archive at the Space Telescope 
Institute (PID-number 8301). STScI is operated by the association of 
Universities for Research in Astronomy, Inc. under the NASA contract 
NAS 5-26555. Also based on observations made with ESO Telescopes 
at the La Silla or Paranal Observatories under programme ID 68.A- 
02555 and 073.A-0050 and on observations with XMM-Newton, an 
ESA Science Mission with instruments and contributions directly 
funded by ESA Member states and the USA (NASA). 



However, the mass determination of galaxy clusters is a very 
difficult task. It is dependent on the method adopted and on 
the validity of the assumptions used to convert observables to 
cluster masses. Currently two methods are widely used: (a) 
From the gas density and temperature profiles measured with 
X-ray observations it is possible to derive an estimate of the 
gravitational mass by assuming spherical symmetry and hy- 



drodynamical equihbrium (see e.g. [Allen et aL 



2001 



lEttori & Lombardil [2003i [Pointecouteau et al.| |2004[ [Pratt &1 



2002 



Amaud 2005 Voigt & Fabian 2006 1. (b) The second method 



is based on gravitational lensing analyses, using either strongly 
deformed background sources (arcs) to constrain the cluster 
mass in the very cluster centre or statistical methods to investi- 
gate systematic shape distortions of background objects to map 
the mass distribution of a cluster (weak lensing method, see 



e.g. Bartelmann & Schneider (2001 1 for a review on this topic). 
The lensing method is affected by the least number of assump- 
tions as it is neither sensitive to the nature of the matter nor its 
dynamical state. However, this method measures the integrated 
mass along the line of sight, which can lead to a bias of too high 
mass estimates (White et al. 2002[ l. Detailed lensing analyses 
were carried out for several galaxy clusters, e.g. CL0024+1654 
( |Kneib et al.|[2003)[Czoske et al.!'"2002), A2218 ( ,Kneib etH 
A1689 dBroadhurst et al. , 2005aibll, A383 (ISmith et al 



1996 



200T) or RX J1347-1145 ('Bradac etaL]|2005]l. 
Unfortunately the mass estimates derived from the differ- 
ent methods can be quite inconsistent. In some clusters 
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there are considerable discrepancies up to a factor of 3, e.g. 
MS0440+0204 (|Uioia et all[T998l l or CL0500-24 ( |Schindler| 
19991 1. |Alienl ( |1998 ') found the lensing and X-ray method to be 



consistent for cooling flow clusters, whereas for non-cooling 
flow clusters the mass discrepancy between the strong lensing 
method and the X-ray based mass determinations can differ by 
a factor of up to 2 - 4. This mainly comes from the fact that the 
inner core of clusters, where strong lensing occurs, is not well 
described by the usual simple models used in X-ray methods, 
which are based on the assumptions mentioned above. The dis- 
crepancies of the weak lensing and the X-ray method seem to 
be much smaller fWu et al.|[T998| l. 

In this paper we present a combined optical. X-ray and 
lensing analysis of Z3146. This cluster of galaxies is located 
at a =10h23m39.6s, 5 = -h04°11'10" (J2000) with a redshift 
of Zci = 0.2906 ( |Schwope et al. 2000 1 and was the subject 
of many previous optical ( e.g. [Crawford et aL [1999; Edge 
etaL] [20^ [Chapman et aTl[2002[ [ Edge & Frayer[]2003tj^Sand 



Ettori et al. 2001 Fabian et al. 2002 Hicks & Mushotzky 



et al.) [ 20()5[) and X-ray investigations (e.g. Edge et al. ]l994| 



2005) 1. This prominent cluster is one of the most X-ray lumi- 



nous systems in the ROSAT Blight Survey ((Schwope et al.| 
20001 hereafter RBS) having an X-ray luminosity of about 



log(Lx) = 45.3erg/s in the 0.5-2 keV ROSAT band. It is part 
of a larger sample of X-ray selected galaxy clusters which is 
described in Sect. |2] and will be given in more detail in a forth- 
coming paper (Kausch et al., in prep.). This paper contains de- 
tailed X-ray and lensing analyses of Z3 146 followed by several 
optical investigations. The X-ray analysis is presented in Sect. 
[3] a description of the optical observations and the data reduc- 
tion procedure used for this investigation is given in Sect. |4] 
Sect. [5] contains a lensing analysis based on weak (Sect. [5.l| 
and strong lensing (Sect. [5^ . Sect.|6] comprises continuative 
optical investigations on the cluster. In section Sect. [7] we sum- 
marize and discuss the results. 

Throughout this paper we use Hq - 70/i7okms"'Mpc"' and 
= 1 - i^A = 0.3. Hence l'.'0=4.36 h'^ kpc for the cluster 
redshift of z = 0.2906. 



2. Description of the ARCRAIDER-Project 

Z3146 is part of the ARCRAIDER sample of galaxy clusters 
(Kausch et al. in prep.). ARCRAIDER stands for ARCstatistics 
with X-RAy luminous meDium rEdhift galaxy clusteRs. The 
project is based on a homogeneous and unique sample of 
galaxy clusters chosen from the RBS ( [Schwope et al. [2000| ), 
a compilation of all X-ray sources with a PSPC count rate 
> 0.2 s As all sources are located at high galactic latitudes 
i\b\ > 30°), the riu values for our clusters are very small 



(hh ^ 7.7 X 10 cm" ). The selected clusters satisfy the 
following criteria: (a) located in the medium redshift range 
0.1 < z < 0.52, (b) an X-ray luminosity > 0.5 x 10"*^ erg/sec 
(0.5-2 keV band), (c) classified as clusters in the ROSAT 
Bright Survey, (d) not a member of the Abell catalogue, and 
(e) visible from La Silla/Paranal (declination 6 < 20°). 
The total sample contains 22 galaxy clusters which were ob- 
served with different telescopes: RBS1316 (RX J1347-1145) is 



[AUen et al.[ [20021 [Gitti & Schiridlerl |2004l Gitti, Piffaretti & 
Schindler 2007, in prep.) and was observed in the U, B, V, 
R and / band with the ESOVLT with the FORSl instrument 
and in the K, band with ISAAC ( [Bradac et aLl |2005]l. AU 
other clusters were observed at least in the V and R band 
either with the SUperb Seeing Imager 2 (SUSI2@ESONTT, 
ESO-filters V#812 and R#813) or with the Wide Field 
Imager (BB#V/89_ES0843 and WFI@ES02.2m, ESO-filters 
BB#Rc/162_ES0844) with usually half the exposure time in 
V than in R. We use the deep R band frame as the primary 
science band for our lensing analysis, whereas the shallow V 
image is used for colour determinations for a rough division 
between foreground and background galaxies. 
As our clusters are the most luminous ones of the RBS, we 
expect these systems to be very massive due to the Lx-M rela- 
tion ( [Reiprich & Bohring er[ [IW9) |S chind ler| (1999)1. Therefore 
it is very likely to find strong gravitational lensing features like 
arcs or arclets in such systems. A similar sample of clusters 



was established by Luppino et al. ( 1999), based on the EMSS. 
In total they found arc(lets) and candidates in ~ 42% of their 
members. As their X-ray luminosity limit was chosen to be 
lower than ours (Lx > 2- 10"*"* erg/sec in the 0.3-3.5 keV regime) 
we expect to detect gravitational arcs in 45-60% of the clusters. 



3. X-ray analysis of Z3146 

3.1. Observation and data preparation 

Z3 146 was observed by XMM-Newton in December 2000 dur- 
ing rev. 182 (PI: Mushotzky) with the MOS and pn detectors in 
Full Frame Mode with THIN filter, for an exposure time of 53.1 
ks for MOS and 46.1 ks for pn. We used the SASv6.0.0 pro- 
cessing tasks emchain and epchain to generate calibrated event 
files from raw data. Throughout this analysis single pixel events 
for the pn data (PATTERN 0) were selected, while for the MOS 
data sets the PATTERNS 0-12 were used. The removal of bright 
pixels and hot columns was done in a conservative way ap- 
plying the expression (FLAG==0). To reject the soft proton 
flares we accumulated the light curve in the [10-12] keV band 
for MOS and [12-14] keV band for pn, where the emission is 
dominated by the particle-induced background, and excluded 
all the intervals of exposure time having a count rate higher 
than a certain threshold value (the chosen threshold values are 
15 counts/100 s for MOS and 20 counts/100 s for pn). The re- 
maining exposure times after cleaning are 52.3 ks for MOSl, 
52.6 ks for MOS2 and 45.7 ks for pn. Starting from the output 
of the SAS detection source task, we made a visual selection on 
a wide energy band MOS & pn image of point sources in the 
FoV. Events from these regions were excluded directly from 
each event list. 

The background estimates were obtained using a blank-sky 
observation consisting of several high-latitude pointings with 
sources removed (Lumb et al. 20021. The blank-sky back- 



the most X-ray luminous cluster known ( Schindler et al. 1997 



ground events were selected using the same selection criteria 
(such as PATTERN, FLAG, etc.), intensity filter (for flare rejec- 
tion) and point source removal used for the observation events; 
this yields final exposure times for the blank fields of 365 ks 
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for MOSl, 350 ks for MOS2 and 294 ks for pn. Since the cos- 
mic ray induced background might slightly change with time, 
we computed the ratio of the total count rates in the high en- 
ergy band ([10-12] keV for MOS and [12-14] keV forpn). The 
obtained normaHzation factors (0.827, 0.820, 0.836 for MOSl, 
MOS2 and pn, respectively) were then used to renormalize the 
blank field data. The blank-sky background files were recast 
in order to have the same sky coordinates as Z3146. For the 
pn data, we generated a list of out-of-time events (hereafter 
OoT) to be treated as an additional background component. 
The effect of OoT in the current observing mode (Full Frame) 
is 6.3%. The OoT event list was processed in a similar way as 
done for the pn observation event file. The background subtrac- 
tion (for spectra and surface brightness profiles) was performed 



as described in Arnaud et al. ( 2002 1. In case of pn the OoT data 
were also subtracted. 

The source and background events were corrected for vi- 



gnetting using the weighted method described in Arnaud et al 



(200111, the weight coefficients being tabulated in the event list 
with the SAS task evigweight. This allows us to use the on-axis 
response matrices and effective areas. 

Unless otherwise stated, the reported eiTors are at 90% confi- 
dence level in the entire Sect. [5] 



We notice a sharp central surface brightness peak at a posi- 
tion 10'>23"40^01 + 04°11'09'.'45 (J2000), in very good agree- 
ment (Aa - 0'.'41, M - 0'.'55) with the optical position of the 
central dominant cluster galaxy (Schwope et al. 2000| l. The 
morphology of the cluster is quite regular, thus indicating a re- 
laxed dynamical state, even though we notice that the central 
core appears slightly shifted to the south-east with respect to 
the outer envelope, with a north-west to south-east elongation 
of the cluster core. The regular morphology of the cluster is in- 
dicative of a relaxed dynamical state, thus allowing us to derive 
a good mass estimate based on the usual assumptions of hydro- 
static equilibrium and spherical symmetry (see Sect. |3.5| l. 

3.2.1 . Surface brightness profile 

We computed a background-subtracted vignetting-corrected ra- 
dial surface brightness profile in the [0.3-2] keV energy band 
for each camera separately. The profiles for the three detectors 
were then added into a single profile and binned such that at 
least a sigma-to-noise ratio of 3 was reached. The cluster emis- 
sion is detected up to 1.5 Mpc (~ 6') and the profile appears rel- 
atively regular and relaxed (see Fig.|2]). The surface brightness 



3.2. Morphological analysis 

The adaptively smoothed, exposure corrected MOS 1 count rate 
image in the [0.3-10] keV energy band is presented in Fig.[T] 
The smoothed image was obtained from the raw image cor- 
rected for the exposure map (that accounts for spatial quan- 
tum efficiency, mirror vignetting and field of view) by running 
the task asmooth set to a desired signal-to-noise ratio of 20. 
Regions exposed with less than 10% of the total exposure were 
not considered. 




Fig.l. MOSl image of ZW3146 in the [0.3-10] keV energy 
band. The image is coiTected for vignetting and exposure and 
is adaptively smoothed (signal-to-noise ratio = 20). 



X 



■Ml 



1?^ 




IP* 



Fig. 2. Background subtracted, azimuthally averaged radial 
surface brightness profile in the [0.3-2] keV range. The best 
fit yS-model fitted over the 230 - 1300 kpc region (0'.9-5'.0) is 
over-plotted as a solid line. When extrapolated to the centre, 
this model shows a strong deficit as compared to the observed 
surface brightness. 



profile of the undisturbed cluster was fitted with the CIAO tool 
Sherpa with various parametric models, which were convolved 
with the XMM-Newton PSF The overall PSF was obtained by 

estimated 



adding the PSF of each camera (Ghizzardi 2001 



at an energy of 1 .5 keV and weighted by the respective clus- 
ter count rate in the 0.3-2 keV energy band. A single /3-model 



4 



W. Kausch et al.: ARCRAIDER I: Detailed optical and X-ray analysis of the cooling flow cluster Z3146 



( [Cavaliere & Fusco-Femiaiiol |1976| l is not a good description 
of the entire profile: a fit to the outer regions (230 /z^q kpc < r < 
1300 h^iQ kpc) shows a strong excess in the centre as compared 
to the model (see Fig. |2|i. The peaked emission is a strong in- 
dication for a cooling flow in this cluster We found that for 
230 /Zyy kpc < r < 1300 /i^g kpc the data can be described 
Of red ~ 1-39 for 74 d.o.f.) by a /3-model with a core radius 
Tc = 177 ± 2 hZi^ kpc and a slope parameter fi - 0.11 ± 0.01 
(3 o" confidence level). The single j6-model functional form is a 
convenient representation of the gas density profile in the outer 
regions, which is used as a tracer for the potential. The param- 
eters of this best fit are thus used in the following to estimate 
the total mass profile in the region where the single beta model 
holds (see Sect. |3.5| ). 

We also considered a double isothermal y6-model and found 
that it can account for the entire profile, if the very inner and 
outer regions are excluded: for 15 kpc h^^ < r < 1300 /ly,,' kpc 
the best fit parameters are = 111 + 2 h:^^ kpc, = 
39 ± 1 h-^ kpc and /? = 0.76 ± 0.01 (xi^ ~ 1.67 for 95 d.o.f.). 
A common /3 value is assumed in this model, but we also tried 
the fit with two different p values, finding very similar results 
(rei = 172 + 2 kpc, r^2 = 45 ± 1 /i^,,' kpc, /3i = 0.76 + 0.01 
and /32 = 0.88 ± O.OUxia ~ 1-68 for 94 d.o.f.). 




Fig. 3. Temperature map obtained by using 4 X-ray colours 
(0.3-1.0, 1.0-2.0, 2.0-4.5, 4.5-8 keV) and estimating the ex- 
pected count rate with XSPEC for a thermal MEKAL model, 
with fixed Galactic absorption Nh and metallicity Z. 



3.3. Temperature map 

The temperature image of the cluster central region shown in 
Fig. [3] was built from X-ray colours. Specifically, we extracted 
mosaiced MOS images in four different energy bands (0.3-1.0 
keV, 1.0-2.0 keV, 2.0-4.5 keV and 4.5-8 keV), subtracted the 
background and divided the resulting images by the exposure 
maps. A temperature in each pixel of the map was obtained 
by fitting values in each pixel of these images with a thermal 
plasma, fixing nu to the Galactic value and the element abun- 
dance to 0.3 solar In particular we note that the very central 
region is cooler than the surrounding medium and the north- 
west quadrant appears slightly hotter than the south-east one, 
even though no strong features are present. The regularity of the 
temperature distribution points to a relaxed dynamical state of 
the cluster, thus excluding the presence of an ongoing merger 
Since cluster merging can cause strong deviations from the as- 
sumption of an equilibrium configuration, this allows us to de- 



rive a good estimate of the cluster mass (see Sect. 3.5 i. 



3.4. Spectral analysis 

Throughout the analysis, a single spectrum was extracted for 
each region of interest and was then regrouped to reach a signif- 
icance level of at least 3 cr in each bin. The data were modelled 
using the XSPEC code, version 11.3.0. Unless otherwise 
stated, the relative normalizations of the MOS and pn spectra 
were left free when fitted simultaneously. We used the follow- 
ing response matrices: ml_169_im_pall_vl . 2 .rmf (MOSl), 
m2_169_im_pall_vl.2.rmf (MOS2), epn_f £2(8_sY9 . rmf 
(pn). 



3.4.1. Global spectrum analysis 

For each instrument, a global spectrum was extracted from all 
events lying within 5' of the cluster emission peak. We tested in 
detail the consistency between the three cameras by fitting sep- 
arately these spectra with an absorbed MEKAL model with the 
redshift fixed at z=0.291 and the absorbing fixed at the galac - 
tic value (hh = 3.01 x 10^° cm"^, [Dickey & Lockman|[T990| . 
Fitting the data from all instruments above 0.3 keV led to in- 
consistent values for the temperature derived with the MOS and 
pn cameras: kT = 6.18+[J |^ keV (MOSl), 5.12+l \l (MOS2), 
5.01;^|j Qg (pn). We then performed a systematic study of the ef- 
fect of imposing various high and low-energy cutoffs, for each 
instrument. Good agreement between the three cameras was 
found in the [1.0-10.0] keV energy range (kT = 6.27+° ?^ keV 
for MOSl, 5.99+!!?[! for MOS2, 6.00+° '^ 



0.20 

Q jg v^.vjvj Q for pn). On the other 

hand, we also found consistent results by fitting the MOS spec- 
tra in the [0.4-10] keV energy range and the pn spectrum in 
the [0.9-10] keV energy range. The discrepancies observed by 
fitting the whole energy range are probably due to some resid- 
ual calibration uncertainties in the low-energy response of all 
instruments. Thus, in order to avoid inaccurate measurements 
due to calibration problems, we adopted the low energy cut-off 
derived above for the spectral analysis discussed below. 
The combined MOS+pn global temperature and luminosity are 
respectively kT = 5.91+[J[;^ keV, Lx (2-10 keV) = 2.0 + 0.1 X 
lO'*^ erg s-i in the [0.4/0.9-10.0] keV energy range (MOS/pn) 



and kT 



6.08tli 



keV, Lx(2-10keV) 



2.0 + 0.1 X 



10"*^ erg s"' in the [1.0-10.0] keV energy range. These values 
are in agreement with ASCA results (Allen et al. 1996, Allen 
2000), while Ettori et al. (2001) derived higher temperature val- 
ues from BeppoSAX observations. The simultaneous fit in the 
[0.4/0.9-10.0] keV energy range (MOS/pn) to the three spectra 
is shown in Fig. |4] 
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Table 1. Results from the spectral fitting in concentric annular regions in the [0.4/0.9-10.0] keV energy range (MOS/pn) and in 
the [1.0-10.0] keV energy range. Temperatures (kT) are in keV, metallicities (Z) in solar units and [2-10] keV luminosities (Lx) 
in units of lO"*"* erg s"'. The total values and numbers of degrees of freedom (dof) in the fits are also listed. Error bars are at 
the 90% confidence levels on a single parameter of interest. 



[0.4/0.9-10.0] keV energy range (MOS/pn) 
Radius Radius (kpc) kT Z Lx ;f^/dof 



0-20" 
20" -40" 
40"- 1' 
1'-1'.5 
l'.5-2' 
2' -3' 
3'-5' 
0'-5' 



0-87 
87-174 
174-262 
262-393 
393-524 
524-785 
785-1309 
0-1309 



4 7+0.1 
^- '-0.1 

6 7+"-3 
6 5+"-3 

U.J 03 

' _o.4 

7.0t« 
6.8!ii 
5.9+!!-! 



0.38!°:°^ 

U./8_oog 
Q 24+0.10 

"•^^-0.10 

23*°-^^ 

"•^-'-0.22 

0.30+!;!!? 



5.68 
5.15 
3.14 
2.41 
1.28 
1.21 
0.79 
20.14 



1630/1503 
1456/1421 
1145/1139 
1002/1007 
743/719 
593/578 
306/253 
2632/1859 



[1.0-10.0] keV energy range 
kT Z Lx x^/dof 



4 8+"' 

"•-'-0.2 

6 9+°^ 

6 7+° '* 

-0.3 

'•"-0.5 

7 4+0.8 

6 r"' 



36+""'' 
30+°"'' 

"•-'"-004 

0.25!-? 
0.18!«;«= 

22+" '° 
"•^ -OlO 

21+" * 

"•^'--021 

0.28+""? 



5.77 
5.20 
3.16 
2.43 
1.30 
1.26 
0.80 
20.42 



1461/1403 
1323/1321 
1040/1039 
892/907 
618/620 
479/480 
241/195 
2399/1759 







■ 




R . ■ 







Fig. 4. Global MOS (lower) and pn (upper) spectra in the 
[0.4/0.9-10.0] keV energy range (MOS/pn) integrated in a cir- 
cular region of radius 5'. The fit with a MEKAL model and the 
residuals are shown. 



3.4.2. Radial temperature profile 

We produced a radial temperature profile by extracting spectra 
in annuli centred on the peak of the X-ray emission. The an- 
nular regions are detailed in Table[T| The data from the three 
cameras have been modelled simultaneously using a simple, 
single-temperature model (MEKAL plasma emission code in 
XSPEC) with the absorbing column density fixed at the nom- 
inal Galactic value. The free parameters in this model are the 
temperature kT, metallicity Z (measured relative to the solar 
values, with the various elements assumed to be present in 
their solar ratios) and normalization (emission measure). We 
separately performed the spectral fitting in the [0.4/0.9-10.0] 
keV energy range (MOS/pn) and in the [1.0-10.0] keV energy 
range. The best-fitting parameter values and 90% confidence 
levels derived from the fits to the annular spectra are summa- 
rized in Table[T] The projected temperature profile determined 
with this model is shown in Fig. |5] We note that, as expected. 



temperature values derived in the [1.0-10.0] keV energy range 
are slightly higher than those derived in the [0.4/0.9-10.0] keV 
energy range (MOS/pn), even though consistent within the 90% 
confidence level. In the following discussion we adopt results 
derived in the [0.4/0.9-10.0] keV energy range (MOS/pn). The 
temperature rises from a mean value of 4.7 + 0.1 keV within 
90 h-^ kpc to kT = 6.7 + 0.5 keV over the 180-1300 kpc region, 
where the cluster can be considered approximately isothermal. 
The lack of evidence for a temperature decline in the outer re- 
gions is in agreement with the results by Mushotzky (2003). 



r [arcmin] 
0.25 0.5 1 2 3 
"1 T 



T" 



T" 



-i i 



[kpc] 



Fig. 5. Full Circles: the projected X-ray gas temperature pro- 
file measured from XMM data in the [0.4/0.9-10.0] keV energy 
range (MOS/pn). Open Triangles: same as full circles, but in 
the [1.0-10.0] keV energy range. 
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The metallicity profile is shown in Fig. |6] a gradient is visible 
towards the central region, the metallicity increasing from Z = 
0.18 ± 0.05 over the 260-400 h'^^ kpc region to Z = 0.38 ± 0.03 
inside the central 90 /i^i kpc. 



r [arcmin] 
0.5 1 



S 0-2 




r [ kpc ] 



Fig. 6. Full Circles: the projected X-ray gas metallicity pro- 
file measured from XMM data in the [0.4/0.9-10.0] keV energy 
range (MOS/pn). Open Triangles: same as full circles, but in 
the [1.0-10.0] keV energy range. 



3.4.3. Cooling core analysis 

The surface brightness profile, the temperature map and the 
temperature profile all give hints of the presence of a cooling 
core. Here we further investigate the physical properties of the 
ICM in the central region. The cooling time is calculated as the 
time taken for the gas to radiate its enthalpy per unit volume 
using the instantaneous cooling rate at any temperature: 



H 



J 



kT 



n^nuA(T) j - 1 ^iX^nMT) 



(1) 



where j - 5/3 is the adiabatic index; // » 0.61 (for a fully- 
ionized plasma) is the molecular weight; Xn » 0.71 is the hy- 
drogen mass fraction; and MT) is the cooling function. We cal- 
culate the cooling function and the electron density by follow- 



ing the procedure described in Sect. 3.5 using the yS-parameters 
derived by fitting the surface brightness profile over the 65 - 
1300 /;^,5 kpc regiorQ(data in this region can be approximated 
by a /? model with r^ ~ 113 /i^^ kpc and a slope parameter 
j3 ~ 0.70). The cooling time is less than 10 Gyr inside a radius 
of 150 /i^y kpc (~ 0'.6), in agreement with the ROSAT result of 
Allen (2000). 

We therefore accumulate the spectrum in the central 0'.6. We 



' The best fit obtained in Sect. 3.2.1 cannot be extrapolated to the 
central region (r < 230 kpc) therefore cannot be used here for the 
purpose of calculating the central cooling time. 



Table 2. The best-fit parameter values and 90% confi- 
dence limits from the spectral analysis in the central 0'.6 
region. Temperatures (kT) are in keV, metallicities (Z) as 
a fraction of the solar value and normalizations in units 
of 10-'4„^„py/4;r[Z)A(i + 2)]2 as done in XSPEC (for the 
MKCFLOW model the normalization is parameterized in terms 
of the mass deposition rate M, in M0 yr"'). 



Parameter 


MEKAL 


MEKAL+MKCFLOW 


kT 
Z 
Norm 

kT\a^ 


5 2+°-' 
6.67:°:»^ X 10-3 


0.40!«;« 

1.07+!J ^^ X 10-3 

1 7+''-2 
^•'-0.2 


Normiow 




M = 1580!;^o 


xVdof 


2024/1783 


1903/1781" 



compare the MEKAL model akeady used in Sect. 3.4.1 and 



3.4.2 with a model which includes a single temperature com- 
ponent plus an isobaric multi-phase component (MEKAL + 
MKCFLOW in XSPEC), where the minimum temperature, 
kTiov,, and the normalization of the multi-phase component, 
Normio^v - M, are additional free parameters. The maximum 
temperature kT of the MKCFLOW model is linked to the am- 
bient value of the MEKAL model. This model differs from the 
standard cooling flow model as the minimum temperature is 
not set to zero. The results, summarized in Table|2] show that 
the statistical improvements obtained by introducing an addi- 
tional emission component compared to the single-temperature 
model are significant at more than the 99% level according to 
the F-test, with the temperature of the hot gas being remark- 
ably higher than that derived in the single-phase model. The 
fit with the modified cooling flow model sets tight constraints 
on the existence of a minimum temperature (~ 1.7 keV). We 
find a very high value of the nominal mass deposition rate in 
this empirical model: ~ 1600 M© yr-'. ASCA-ROSAT obser- 
vations already found a very strong cooling flow in this cluster 
(Allen 2000). 

3.5. Mass determination 

In the following we estimate the total mass of the cluster using 
the usual assumptions of hydrostatic equilibrium and spheri- 
cal symmetry. Under these assumptions, the gravitational mass 
Mtot of a galaxy cluster can be written as: 



Mtot(< r) = - 



kTr 
Gfinip 



d\nn„ dlnT 

+ 

dlnr dlnr 



(2) 



where G and nip are the gravitational constant and proton mass 
and fi X 0.61. The deprojected dlnng/dlnr was calculated 
from the parameters of the single j6-model derived in Sect. 



3.2. 1 In particular, the advantage of using a y6-model to param- 
eterize the observed surface brightness is that gas density and 
total mass profiles can be recovered analytically and expressed 
by simple formulae: 



ng3s(r) = no,! 



1 + 



-3/3/2 



(3) 
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M,o,{< r) = 



Gfinip 



3y6rr dT 



dr 



(4) 



In estimating the temperature gradienj^from the profile shown 
in Fig. [5] only data beyond 30" (~ 130 h^l kpc) were con- 
sidered: in the central bins the temperature as derived in Sect. 



3.4.2|is more aff'ected by the XMM PSF and projection effects. 



while for the outer regions these effects can be neglected (e.g. 
Kaastra et al."2004). The total gravitating mass distribution de- 



rived from Eq. |4]is shown in Fig.|7]as a solid line, with errors 
coming from the temperature measurement and j6-model repre- 
sented as dashed lines. Within 1 Mpc we find a total mass of 
~ (5.4+0.5) X lO'^^Mo and within the outer radius of the cluster 
as visible in the X-ray surface brightness profile (1.5 Mpc) we 
find ~ (8.5 ± 0.7) x W'^Mq. We note that the total integrated 
mass within a particular volume is dependent upon the local 
physical properties (temperature and density gradients) and is 
unaffected by the regions interior, or exterior, to that radius. 
The mass profile derived with this method is thus reliable in 
the region where the /3-model is a good representation of the 
observed surface brightness profile (230 kpc < r < 1300 kpc, 
see Sect. 2.3.1), whereas it cannot be extrapolated to the central 
region. 

We also calculate the projected mass along the line-of-sight 
within a cylinder of projected radius r. The integration was per- 
formed out to a radius of ~ 5 Mpc from the cluster centre. The 
projected total mass is shown in Fig.[7]as a dashed-dotted line. 
In Fig. |7] we also show (dotted line) the gas mass profile de- 
rived by integrating the gas density given by Eq. [3]in spherical 
shells and using the yS-model parameters determined in Sect. 
|3.2.1| The normalization of Eq. |3] is obtained from the com- 
bination of the best-fit results from the spectral and imaging 
analyses, which allows us to determine the conversion count 
rate - flux used to derive the bremsstrahlung emissivity that is 
then integrated along the line-of-sight and compared with the 
central surface brightness value. We note that, since we adopt 
the parameters of the yS-model fit in the outer regions, the de- 
rived central electron density («() e ~ 1.7 x 10"^ cm"^) is that 
predicted by the extrapolation of the j6-model fit to the centre 
(see Fig. |2]l. This procedure is nonetheless reliable in estimat- 
ing the gas mass for r > 230 /z^,,' kpc, shown in Fig.|7] 
In order to allow a direct comparison with our weak lensing 



studies in Sect. 5. 1 and to derive an estimate of r2oo we perform 



a fit to the NEW profile ( [Navarro et aT] [19%! [1997] ) given by 



Mdm(< r) = ^nrlp,- \ [ 

3 ln(l + c) - c/(l + c) 



(5) 



where pc is the critical density. The scale radius and the 
concentration parameter c are the free parameters. The best 
fit parameters that minimize the of the comparison be- 
tween the mass predicted by the integrated NEW dark mat- 
ter profile and the mass profile reconstructed from Eq. |4]are: 
Tj = 516.6 + 13 kpc and c - 3.58 + 0.07, where the relation 
'"200 - cx Ts holds. The quoted error are at the 68% confidence 



^ As a first-order approximation, the temperature gradient is esti- 
mated by least-squares fitting a straight line to the observed depro- 
jected temperature profile. 




Fig. 7. Solid line: Integrated total mass calculated from Eq. [4] 
(with eiTor on the mass calculation coming from the temper- 
ature measurement and y6-model shown by the dashed lines). 
Dashed-dotted line : Projected total mass with error Dotted 
line: Gas mass. 



levels (Icr) on a single parameter of interest. Note that we ne- 
glect the gas mass contribution to the total mass and we as- 
sume Mtot(< r) = Mdm(< '")■ We have also performed the same 
fitting procedure by including the gas mass, i.e. by assuming 
Mtoti< r) - Mdm(< f) + Mgas(< r), and find very similar re- 
sults. In Fig.[27|a) we show confidence contours of the NEW 
fit to our mass profile and note that both parameters are well 
constrained. For a comparison to the equivalent model based 
on lensing data we refer to Sect.|7] 

4. Optical observations 

4.1. WFI-observations and data reduction 

Z3146 was observed with WFI@ESO/MPG2.2m in the two 
observing programs 68.A-0255 (P.I. S. Schindler) and 073.A- 
0050 (PI. P Schneider). The first one obtained 8000s in broad 
band V (BB#V/89_ES0843) and 16100s in broad band R (ESQ 
filter BB#Rc/162_ES0844). The second programme observed 
for another 8900s in R and 1500s in B (BB#B/123_ES0878). 
All data were processed with our image reduction pipeline de- 
veloped within the GaBoDS project. It performs all necessary 
steps from raw frames to astrometrically and photometrically 
calibrated and co-added images. The individual methods and its 



performances on WFI data are described in detail in Schirmer 



let al.| P003| l and |Erben et alT] ( |2005| . AU data were obtained 
during clear nights, under good seeing conditions (< 1'.'2) 
and with a large dither box of about 3'.0 to ensure good flat- 
fielding and an accurate astrometric calibration. The images 
were tied to the astrometric frame of the USNO-A2 catalogue 



8 



W. Kausch et al.: ARCRAIDER I: Detailed optical and X-ray analysis of the cooling flow cluster Z3146 




- #1 



, I 
- 1 



Fig. 8. Layout of the Z3146 R-hmd WFI observations of the 
Bonn (B; solid lines) and Innsbruck (I; dashed lines) groups. 
The orientation with respect to the sky is North-Up and East- 
Left. The overlap in the R exposures allows useful tests in the 
shape measurement procedures during our weak lensing anal- 
ysis (see Sect. |5.1[ ). If not stated otherwise we use the deep 
R band stack from both observing campaigns and work in the 
common area (A) where we have data in B, V and R available. It 
has an extent of 25'.8 x 23'.8 and the cluster is located 8'.0 north 
and ICXO west with respect to its lower left corner. The white 
square marks the ~ 16'40" x 16'40" region used in Sect. 5.2 
and[6]for extracting the Red Sequence galaxies. 



(Monet et al. 1998[), pho tometrically calibrated with Stetson 
standards (Stetson 2000) and finally co-added with the Swarp 
tool ( |Bertin,,2002) . We produced several co-added images from 
our R band exposures mainly to crosscheck our object shape 
measurements in the weak lensing analysis (see Sect. |5.1| l. The 
characteristics of all co-added images used in this work are 
summarised in Table[3]and Figs. [8] and |9] Each co-added sci- 
ence image has a pixel scale of 0'.'238 and is accompanied by a 
weight map characterising its noise properties. 

4.2. HST-observations 

In addition to our WFI observations we used archival calibrated 
HST data from the WFPC2 Associations obtained during a 
snapshot programme (PI: Edge, PID-number 8301). Z3146 
was observed with the WFPC2 (filter F606W) in April 2000 
for a total exposure time of 1000 s. The pixel scale is ~ (X'l per 
pixel, the seeing is measured with the FWHM_IMGE keyword 
of SExtractor to be ~ 0.1". The main purpose using these 
archival HST data was the identification and comparison of arc 
candidates with the ground based observations. 




Fig. 9. As a control check for our photometric calibration we 
plot colours of bright, unsaturated stars in our WFI fields com- 
pared to isochrones for stars of 10 Gyrs with a metallicity 
of 0.4% (see [Girardi et aL) |2002| ). Magnitudes for the stars 
were estimated with the SExtracor MAG^UTO magnitude. 
Our data are in excellent agreement with the isochrone predic- 
tions. 



5. Lensing analysis 

5.1. Weak lensing analysis 

As a second method for determining the mass and its distribu- 
tion of Z3146 we performed a weak lensing mass reconstruc- 
tion. For a broad introduction on weak lensing and its applica- 
tion in cluster mass determinations see for instance iBartelmarml 



& Schneider (2001 1. In the following we describe the creation 



of our background galaxy catalogue and the weak lensing anal- 
ysis. Throughout the analysis we use standard weak lensing 
notation. 

5.1 .1 . Lensing catalogue generation 

We use the deepest image, the /?-band for the weak lens- 
ing measurement. To create an object catalogue with shear 
estimates for all objects we first extract sources with the 
SExtractor (we use a detection threshold of 1.9 and a min- 
imum area of 3 pixels for our detections) and the IMC Al]^ soft- 
wares. While SExtractor produces a very clean object cata- 
logue if the source extraction from the science images is sup- 



ported by a weight map (see e.g. Fig. 27 of Erben et al. 2005 1, 
IMCAT calculates all quantities necessary to estimate object 
shapes. We merge the two catalogues and calculate shear esti- 



mates with the KSB algorithm (see Kaiser et al. 1995 i. For the 
exact application of the KSB formalism we closely follow the 
procedures given in Erben etaL] |2001^ with important mod- 
ifications in the selection of stars that are used for the neces- 



see http://www.ifa.hawaii.edu/~kaiser/imcat 
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Filter 


Image Region 


exp. time (s) 


limiting magnitude (3ct in a 2'.'0 aperture) 


Seeing (arcsec) 


B 


(B) 


1500 


24.70 


1.21 


V 


(I) 


8040 


25.02 


1.16 


R 


(I) 


16077 


25.27 


1.02 


R 


(B) 


8850 


25.24 


1.11 


R 


(A) 


24927 


25.71 


1.04 



Table 3. Characteristics of the co-added WFI images. The limiting magnitudes quoted in column 4 are defined in the Vega 
system via nin^ - ZP - 2.5 log( ^JW^ ■ 3 ■ cr), where ZP is the magnitude zeropoint, A^pix is the number of image pixels in a circle 
with radius 2'.'0 and cr the sky background noise. The seeing in column 5 was measured with the SExtractor FWHM .IMAGE 
parameter. The groups in Bonn (B) and Innsbruck (I) observed Z3146 with an offset of about 5'.0 In Ra and 8'.0 in Dec. The exact 
layout is given in Fig. [8] For the R band we created, besides the individual stacks, a deep mosaic in the common area (A). We 
note that the Innsbruck R band has nearly double the exposure time of the Bonn image but about the same limiting magnitude as 
it was observed during less favourable moon phases. The quality of our photometric calibration is crosschecked in Fig. [9] 



sary PSF corrections (see Van Waerbeke et al. 2005[ l. After 
the PSF corrections we reject all objects with an IMCAT signifi- 
cance y < 8, a half light radius smaller or equal to that of bright 
stars, and a final modulus of the shear estimate \e\ > 0.8. In the 
following analysis we do not apply a weighting to individual 
galaxies. 

Because of the large offsets between the observations from 
the Bonn and Innsbruck groups it is not obvious whether we 
can savely use the deep stack (A) for the lensing analysis or 
whether we have to work on the individual co-additions (I) and 
(B). This mainly comes from the assumption of a smooth vari- 
ation of the PSF over the whole field-of-view when we correct 
galaxy shapes for PSF effects. With the large offsets we could 
suffer from discrete jumps in the PSF anisotropy within the (A) 
mosaic. However, in Fig. 10 we see that the PSFs of all three 



stacks are well behaved in the (A) area. We performed compar- 
isons of the final shear estimates in the (I), (B) and (A) mosaics 
and found that they are in excellent agreement (see Figs.[TT|and 
\T2\ . Therefore we decided to use the (A) mosaic, which is the 
deepest image, for our primary analysis. 

The next step is to clean our lensing catalogue from likely 
cluster members, foreground galaxies and faint stars. To this 
end we plot stars and bright galaxies (R < 22.0), which have a 
high probability to be at a lower redshift than the cluster, in a 



colour-colour diagram (see Fig. 13 i. We note that all these ob- 
jects mostly occupy a limited and well defined area in colour- 
colour space. We finally use the following criteria to clean our 
background source catalog: We reject all objects with R < 22.0 
and keep objects between 22.0 < R < 23.0 if they do not 
lie in the following area of the (B - V) vs. (V - R) diagram: 
-0.23 < (V - R)- 0.8(B -V) < +0.8; 0.25 < (V - R) < 
1.67; 0.2<(B-y)< 1.5. We keep all objects with 7? > 23.0 as 
probable background sources. We note that |Clowe & Schneider] 
( [200 2') and'Dietrich et al. ( 2005 1 used similar criteria to identify 
foreground objects. 

For our final lensing catalogue we additionally exclude all 
objects falling in masked image regions (around bright stars, 
satellite tracks etc.). This leaves us with about 12 galaxies per 
sq. arcmin as direct tracers of the cluster shear. Around the 
Brightest Cluster Galaxy (BCG) of Z3146, this average den- 




Fig. 10. Stellar anisotropics for the fields (I), (B) and (A) in 
their common area.The length of the bars gives a measure for 
the stellar ellipticity and the direction the orientation of the 
semi-major axis. The longest bars coiTespond to an ellipticity 
of about 5%. All patterns show a smooth variation over the 
whole field-of-view of the (A) area. 



sity is reached at a radius of about 45'.' with only very few 
sources in the cluster centre. 



5.1 .2. Weak lensing mass determination of Z3146 

The main interest of the weak lensing analysis in this work is 
an estimate of the total cluster mass of Z3146 and its inter- 
comparison with the X-ray analysis. We first perform a stan- 



dard KS93 cluster mass reconstruction (see Kaiser & Squires 



|1993| ) to investiagte the dark matter distribution and to obtain 
an estimate for the cluster centre. In addition, we calculate a 
B-mode map, i. e. we performed another mass reconstruction 
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Fig. 11. Galaxy ellipticities after all PSF corrections from com- 
pletely independent analyses of the (I), (B) and (A) images in 
the overlap area. In the upper panels we compare the (B) and (I) 
sets which are independent (in terms of raw images entering the 
co-added stacks) from each other. The slopes of the fitted lines 
are 0.983 for el and 0.95 for e2. For the lower panels (compar- 
ison between (A) and (1)) the slopes are 0.992 for el and 0.974 
for e2. The corresponding numbers for (A)-(B) are 1 .080 for e 1 
and 1.091 for e2 (not shown in the figures). We conclude that 
the measurements are in very good agreement with each other. 
For the fits, the Ordinary Least Squares bisector method (see 
|Isobe et aL]|1990[ ), which considers errors in el and e2, is cho- 
sen. After applying the line fits to the ellipticity components 
the mean shear differences are compatible with zero in all four 
cases. We note that the standard deviation of the residuals cr^e 
is 0.16 in both components for the (B)-(l) comparison but 0. 10 
in the (A)-(I) case. The cause for the lower value in the (A)-(I) 
case is that the data are not independent but share a large frac- 
tion of the input images. The residuals as function of position 



are discussed in Fig. 12 



after all object ellipticities have been rotated by 45 degrees. 
This map should contain noise only if the lensing data are free 



from systematics. The results are discussed in Fig. 14 We see 
that our lensing centre is in excellent agreement with that de- 



termined from our X-ray analysis (see also Sect. 3.2 1 and we 
use the X-ray position for the following analysis. We estimate 
significances for peaks in our reconstructions and errors on the 
lensing centre with the following procedure: We randomise the 
orientation of our galaxies, redo a KS93 mass reconstruction 
with the new catalogue and repeat this procedure many times. 
For the peak significance we count how often the k value in our 
noise maps exceeds that of the lensing signal. With 29700 real- 
isations the probabilities that the cluster peak, the cluster mass 
extension and the eastern and western holes in the B-mode map 
are pure noise features are 0/29700, 170/29700, 88/29700 and 



Fig. 12. Differences of coiTected galaxy ellipticity measure- 
ments from the (1) and (A) images as function of image po- 
sition. We see no systematic trends of residuals with object po- 
sition. Also a closer, visual inspection of the two-dimensional 
ellipticity residual distribution reveals no systematics on small 
scales. The standard deviation of the ellipticity components 
around zero is cr^e ^ 0.1 in all four cases. Hence we conclude. 



together with the results of Fig. 1 1 that it is safe to use the (A) 
image for our weak lensing analysis. 



9/29700. Assuming Gaussian statistics this translates to signif- 
icances of > 4.1cr, 2.7cr, 2.9cr and 3.5cr respectively. We con- 
clude that the highly significant cluster peak has no significant 
extension to the South. Next we estimate errors on the lensing 
peak position. The best way to measure the centroid dispersion 
would be to use a parametric model for the mass concentration 
and to generate noisy mass maps with randomised ellipticity 
orientations and galaxy positions. If the model were true we 
would obtain accurate eiTor estimates from the noisy mass real- 
isations. With the data at hand we can follow this idea by con- 
sidering the original reconstruction as the input mass model. 
We probably overestimate the true error in this way because 
our input model already contains measurement noise. We plot 
the result of this exercise for 200 maps in the lower left panel 



of Fig. 15 We find a significant asymmetry in the distribution 
of positional differences. The positional accuracy is about 2.3 
times better along R.A. than in Dec. As the mass distribution 
is elongated towards the South we would expect a skewed dis- 
tribution of the positional errors towards negative Dec values 
but the observed symmetric elongation in the North-South di- 
rection is surprising. We checked that not a few, very elongated 
galaxies or shot noise from the galaxy positions (introduced for 
instance by object masks) are responsible for this result (see 
Fig. [TSj l. Given this result we quote the positional accuracy of 
our lensing centre as ARa - 19'.' 8; ADec - 45'.' 6. 
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Fig. 13. To obtain criteria to reject probable low z foreground 
objects from our lensing analysis we plot all bright objects with 
R < 22.0 in a colour-colour diagram. Most of them lie in a well 
defined area which we use to clean our catalogue. Elliptical 
cluster galaxies for a z = 0.3 cluster are expected around B - 
V ~ 1.2; V-R~ 1.2. 



While our mass maps give us insight into the dark mat- 
ter distribution in Z3146 it is difficult to obtain reliable esti- 
mates for the total lensing mass and the involved errors. The 
main problems are that mass reconstructions involve a convo- 
lution from the measurable shear field and, in addition, they 
become very noisy at a distance several arcminutes from the 
cluster centre (see Fig. 



14 1. Moreover, they intrinsically suf- 



fer from the mass-sheet degeneracy. Hence, we will estimate a 
lensing mass by directly fitting parametrised lensing models to 
the shear data. The error analysis is simplified significantly in 
this case. Moreover, model fits break the mass-sheet degener- 
acy by the explicit assumption that k at large distances from the 
lensing centre is zero. The main drawback is that shear data 
alone do not allow a clear discrimination between different. 



plausible mass models (see Fig. 16 1. 

For our model fits to the shear data we primarily consider 
the universal density profile (NFW) proposed by|Navarro et al.| 
1996| l. The details for the calculation of the lensing quantities 



K and 7 for this profile are given in several publications and 
the details are not repeated here (see e.g. |Bartelmann] |1996[ 
Kruse & Schneider 1999 1. To determine our model parameters 



we use the log-likelihood method proposed in Schneider et al. 
( |2000| l. We maximise the likelihood function: 



exp 



,fl)P 



(6) 



the model, a the set of parameters to be fitted, 6 the galaxy po- 
sition in the lens plane and cr the dispersion of the observed 
ellipticity. It is given by cr = (1 - \g\^)cr^^, where cTe^ stands 
for the (two-dimensional) dispersion of the unlensed source el- 
lipticities. All the model parameters a are contained in g. The 
fundamental assumption of this method is that the source ellip- 
ticity distribution can well be described by a Gaussian distri- 
bution of width cr^^. It is optimal in the sense that it uses the 
full ellipticity information and not only individual components 
(such as fits to the tangential part of the shear). For a more de- 
tailed discussion on this likelihood method and its properties 
see |Schneider et al.| ( |2000l i. 

Before we apply this method we still have to specify the 
redshifts of the source galaxies and the galaxy sample we in- 



clude in our fits. In Fig. 16 we show the tangential component 
of the shear around the cluster centre. We can trace the cluster 
shear signal over the whole field-of-view of 15'.(Q Hence, we 
include all preselected background galaxies in our estimations. 
For the dispersion of the unlensed ellipticities cr^^ we use the 
measured value of our galaxies cr^, averaged over the whole 
field-of-view. Here we assume that weak lensing does not 
change this value significantly and we estimate o-,,^ = 0.38. For 
the redshift distribution of our background galaxies we use es- 
timates from [Hetterscheidt et al.| ( [2006| l. The authors obtained 
photometric redshifts for 62000 galaxies with 21.5 < ^ < 24.5. 
Their WFI data consist of 1.75 sq. degree of deep UBVRI 



photometry in three different patches (see Hetterscheidt et al^ 
|5006 ; Hild ebrandt etaL] [20061 for details on the data). The pho- 
tometric redshift distribution is parametrised by the following 

j2006l ): 



function introduced by Brodwin et al 



pmiz) = A [piiz) H(zt -z)+ P2(z) H(z - zd] , 
where H denotes the Heaviside step function. 



Pi(z) = 
and 

P2(z) = exp 



exp 



pdzt); 



the normalisation A is obtained by 
dzpliz) + J dzpiiz) = 1. 



(7) 



(8) 



(9) 



(10) 



For Zt = 1 the reported fit parameters are: zo - 0.27, a - 
2.84, /3 = 1.40, y = 2.34 and zi = 2.16. With this choice the 
mean redshift for galaxies behind our cluster at Zd = 0.29 is 
(zbg) - 0.79. We use this distribution to estimate the average 
geometrical lensing factor: 



dzpfit(z) ^ , , , (11) 

-0.29 D,{Z) 



where D^, D^^ and represent angular diameter distances 
observer-cluster, cluster-background source and observer- 
background source respectively. 



where is the number of galaxies, e the observed, two- * We note that for radii larger than 8'.0 we cover the area around the 
dimensional elUpticity of each galaxy, g the reduced shear of cluster completely only North- West of its centre (see Fig.[8| 
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For our fits to the NFW profile, we consider the concen- 
tration c and the radius ^200 (see [Navarro et al.| |1996[ ) as free 
model parameters. With our setup, the application of our pre- 
scription to the shear data leads to best fit values of r2oo = 
1661^325 /j^o'kpc and c - 3.6^2 4- The model has a significance 
of 4.35 over one with zero mass and the errors on r2oo and c 
are at the 90% confidence level. They were estimated with our 
likelihood analysis by keeping c or r2oo at its best fit value and 



leaving r2(K) or c as the only free parameter. In Fig. 27 b) we 
show confidence contours of our analysis and note that both 
parameters are reasonably well constrained except for low val- 
ues of c see Sect.|7]for a comparison to an NFW model based 
on X-ray data). In addition to the NFW profile we also mod- 
eled our shear data by a Singular Isothermal Sphere (SIS) char- 
acterised by its velocity dispersion cr,,. Our best fit model has 
cr,, - 869^j53km/s. The errors represent the 90% confidence 
level and the model has a significance of 4.88cr compared to the 
zero mass model. In contrast to the NFW fit, the significance 
and the estimated velocity dispersion of the SIS model show 
some dependence on the galaxies included close to the cluster 
centre. We notice an increase of significance by excluding the 
galaxies in a circle 30'.'0 around the cluster (five objects) and a 
smooth decrease of S /N if we reject more galaxies beyond that 
point. Hence, we used all galaxies with a distance greater than 
30'.'0 from the cluster centre for our SIS fit. 

We finally discuss a possible bias of our result due to 
a systematic underestimate of the shear. As we showed in 
Erben et al.|(|2001) and within the Shear Testing Program (see 



Heymans et al. 2005) our pipeline may underestimate weak 
shear by 10%-15%. We recalculated the best fit NFW values 
after boosting all ellipticities by a factor of 1.15. We then ob- 
tain ^200 - 1'701;^3Q3 /'^flkpc and c - 3.66;^2 ig which is well 
within the error bars of the original signal. Hence, a possible 
systematic underestimate of the shear by about 15% would not 
change our results significantly. At the end of this section we 



show in Fig. 18 the total mass properties given by our model 



fits. We also present a mass-to-light ratio analysis in Sect. 6.3 
and will compare our results with masses from X-ray analyses 
in Sect|7] 

5.2. Strong lensing analysis 

5.2.1. Definition/identification of gravitational arc 
candidates 

In addition to the weak lensing analysis we have searched the 
central cluster region for strongly lensed objects. In ground 
based observations usually only arcs tangentially aligned with 
respect to the mass centre are visible, as radial arcs are very 
thin and faint structures in the vicinity of bright central galax- 
ies of clusters. In addition, arcs and their counter images have 
the same spectra and redshifts > 2 x ziens- However we do not 
have spectra, hence apart from the position and the morphol- 
ogy, the redshift is the main identification criterion. Therefore 
we investigated whether it is possible using our observations to 
roughly estimate the photometric redshift or at least to find out 
whether an object belongs to a fore- or background population. 
For that purpose we have performed simulations using the soft- 
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Fig. 15. The lower left plot shows the lensing centres of our 
cluster reconstitictions after adding noise realisations to our 
signal (see text for details). The quoted distances are with re- 
spect to the original cluster centre. Formal Icr positional un- 
certainties of the lensing peak derived from this distribution 
are AR.A. - 19'.'8 and ADec - 45'.'6. We investigated whether 
the observed asymmetry of the error distribution comes from 
shot noise in the galaxy ellipticities (the error might be domi- 
nated by a few very elliptical galaxies close to the cluster cen- 
tre) or from the galaxy positions (as areas around bright stars or 
other image defects have been masked out the object distribu- 
tion in our field is not homogeneous). To this end we assigned 
each galaxy a random ellipticity of a Gaussian distribution with 
cTj. - 0.38 (see Sect. 15.1.2] ) and repeated our eiTor analysis with 
200 new noise realisations from this catalogue (upper left plot). 
For the upper right plot we additionally randomised the posi- 
tions of our objects. Both simulations show a similar asymme- 
try as the original analysis and we conclude that it does not 
originate from shot noise of the galaxies. We note that the po- 
sitional error estimates of this analysis are an upper limit as the 
original signal is already noisy (we implicitly assumed it to be 
noise free in our calculations). 



ware package hyper ^ ( jBolzonella et al. 20001. We created 
a set of 3000 artificial galaxies with the following parameters: 
< Zsim < 2, R magnitude 22 < niR < 25 (which coiTesponds 
to the range of the arc candidates), using the simulated filter 
WFI R band (ES0844) as the reference filter. The type of the 
galaxies was also randomly chosen to be either E, SO, Sa, Sb, 
Sc or Sd. The simulations have shown that it is not possible 
to obtain any reliable redshift estimate from BVR images only. 
38% of all simulated galaxies with z < Zz3i46 were found to be 
background objects. On the other hand, 27% of the background 
galaxies (defined as z > Zz3i46) were measured to be located in 
front of Z3146. Hence it is even not possible to decide whether 
an object of unknown redshift is a foreground or a background 



http : //webast . ast . obs-mip . fr/hyperz 
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Fig. 14. KS93 weak lensing cluster mass reconstruction of Z3146 (left panel). The shear field was smoothed with a Gaussian of 
width 1'.4. The k map was normalised so that the mean of k in the entire field is zero. Because we use this map only to define the 
centre of mass for model fits to the shear data the actual choice for this normalisation is not important. Solid contours increase 
in steps of Ak = 0.01 (~ 4.29 x lO'^Mg/Mpc^), assuming a source galaxy redshift of 0.79. The cluster clearly shows up as the 
highest peak in this reconstruction. Dashed contours represent the same levels of negative k. The white cross marks the cluster 



centre of our X-ray analysis (see Sect. 3.2 1 which is in excellent agreement with the lensing centre. The peak at the cluster 



position is detected with a significance of more than 4.1cr, the extension to the South (marked with a white circle) with 2.7cr. 
Hence, this extension of the k map is not considered a significant lensing feature. The right panel shows a B-mode map of the 
K field. It contains two deeper holes in the southern part that are significant with 2.9cr (eastern hole) and 3.5cr (western hole). 
Around the cluster position this map shows a null signal. See the text for more details. 
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Fig. 16. Left: Tangential shear signal as function of radius 
from the centre of Z3146. The signal is robust and we can 
trace the cluster shear up to the border of our data field. The 
solid line shows the reduced shear of our best fit NFW model 
('"200 - 1661 /j^gkpc, c - 3.6), the dashed line that of an SIS 
fit (cr,, - 869km/s). See the text for further details. Right: We 
show the cross component of the shear gr around the cluster 
centre, i.e. the signal after rotating all galaxies by 45 degrees. It 
should vanish if lensing caused the original signal and our mea- 
surement is compatible with zero over the entire distance range 



(see also Fig. 14 1. As a further consistency check we compare 



in Fig. 17 tangential shear measurements from the (I), (B) and 



(A) data sets. 



object and we have to restrict our search for strongly lensed ob- 
jects to morphological criteria only. 



Unfortunately there is no common definition of an arc candi- 
date. The definition we adopt of a gravitational arc candidate 
is that of an elongated object, aligned tangentially with respect 
to the cluster center, a minimum length of I'.'O and a length- 
to-width ratio l/w > 1.5. However it is not yet clear whether 
Z3146 can produce strong lensing or not: the low concentra- 
tion parameters c obtained during the modelling of both, the 
X-ray (c = 3.58 + 0.07, rzoo = 1849 ± 47 /17'kpc, see Sect.llsll 



and the weak lensing data (c - 3.6;^2i, r2oo = 1661;^32o h^lkpc, 



see Sect. 5. 1 .2 1 to an NFW profile, leads to an Einstein ring of 
only ~ 1". However, due to the large errors in both the c and 
r2oo determination and the unknown source redshift, we cannot 
exclude the strong lensing ability: a source redshift of z = 2 and 
adopting the upper limits of c and r2oo (leading to c = 6.4 and 
'"200 = 1941 /i^ykpc, based on weak lensing values) shifts the 
Einstein ring to ~ 13". Additionally, our adopted SIS model 
derived in Sect. 15. 1.2] leads to an Einstein radius of the same 
size, assuming the source located at the derived mean redshift 
(Zsource) - 0.79. Adopting cr,, - 993 km/s, the upper limit, leads 
to a critical curve at ~ 22" . Hence we restrict our search to re- 
gions within a radius of about 30", centred on the position of 
the Bright Central Galaxy. 



In a deep arc search using the WFPC2 archive. Sand et al. 
( 2005[ l quote one arc in this archival HST data set (Al in our 
data set). Our identification of strong lensing features was done 
by visual inspection of the WFPC2 frames in direct compari- 
son with the deep WFI exposures. As some of the candidates 
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Fig. 19. HST image of Z3146 obtained by the WFPC2. The exposure time was 1000 s using the filter F606W. The objects A3 



and A4 were found on the HST frame by visual inspection (see Sect. 5.2. 1 for more details). 



are very similar to not fully removed cosmics we carefully 
searched for all identified objects whether there is a corre- 
sponding object on all WFI frames. In this way we identified 4 
objects in total (denoted as Al , . . . , A4) in the chosen field as 



good candidates for being strong lensing features (see Figs. 19 
andra 



A comparison even with shallow space based observations is a 
good method to identify possible gravitational arcs due to the 
missing atmospheric blurring effects. Several of the arc can- 
didates were smoothed on the WFI images so as to even lose 



their tangential alignment. In particular, objects A3 and A4 are 
so strongly influenced by observational effects that they are not 
identify able as arcs in ground based observations. A detailed 
comparison between the WFI and WFPC2 images of the arcs 



is shown in Fig. 20 Note that the exposure time of the WFI R- 
image is 6.9 h, whereas for the WFPC2 it was only ~0.28h. 
Nevertheless, the arc candidates visible in the HST image are 
clearly recognizable as possible gravitational arcs, whereas in 
the WFI frame seeing effects dominate the shape of the objects. 
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Fig. 20. Detailed comparison of the arc candidates in Z3146 (top: WFI, bottom: WFPC2). In all images North is up. East to the 
left, the FoV is 15" x 15". See text for more details. 
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Fig. 17. The figure shows comparisons of the tangential shear 
signal around Z3146 from common background objects for the 

(A) and (B) (upper panel) and the (A) and (1) sets (lower panel). 
The left measurements are from (A) and the shifted ones from 

(B) and (1) respectively. The catalogues of common sources 
were created by merging our lensing objects from (A) with the 
ellipticity catalogues from (1) and (B); see Fig. 11 As in the 



ellipticity comparisons, the tangential shear signals around the 
cluster agree reasonably well. A very different behaviour is ob- 
served for the first bin to which only about a dozen galaxies 
contribute. Note that the initial catalogues were created inde- 
pendently and hence the object samples in the two comparisons 
are not exactly the same. 



5.2.2. Determination of the length-to-width ratio 

To measure the length-to-width ratio l/wwe used SExtr actor 
to detect the arc candidates on the WFPC2 image as it is 
not affected by atmospheric blurring. Due to its shallow- 
ness we used a value of 0.75 for DETECT_THRESH and 
ANALYSIS _THRESH. The l/w ratio itself was determined us- 



ing the same ansatz as in Lenzen et al. ( |2004) and Bertin 
( 2005| l: we treat the arcs as a set of pixels with a certain light 
intensity value at each pixel. The light distribution of a certain 
object is then defined by all corresponding pixels detected by 
SExtractor shown in the SEGMENTATION images. Hence 
we can compute the second moments Ai and A2 of this light dis- 
tribution in the usual way (see e.g. |Lenzen et al. |2004, ,Bertin[ 
|2005[ ). Although the length I is not equal to Ai and the width w 
is not equal to A2 the ratio l/w is equal to A1/A2 ( |Jahne[[2002[ ). 
Hence we obtain the length-to-width ratio by determining Ai 
and /I2. 

5.2.3. Photometry / catalogue creation 

The photometry was also performed with the software pack- 
age SExtractor 2.3.2. In contrast to the determination of 
l/w we used the WFI frames for this purpose, as those images 
are much deeper (see Table|3]l. The photometric measurement 
on the WFPC2 image was skipped as the F606W filter is fully 
covered by the V and R band of the WFI observations. 
As we concentrate on the cluster itself we restricted the ex- 
traction of object catalogues to a FoV of ~ 16'40" x 16'40" 
(4.26 Mpc X 4.26 h^Q Mpc in our cosmology). The V and 
R images were convolved with a slight Gaussian filter of width 
0.61 and 0.91 pixels, respectively, to bring all observations to 
the same seeing of ~ l'.'2 (AFWHM < 10"^") and hence en- 
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Fig. 18. Total weak lensing mass calculated from our best fit 
mass models in spheres of radius r around the cluster cen- 
tre. Solid lines encompass our best fit NFW model r20() = 
1661^325 /Zyglcpc; c - 3.6 with the errors from the mea- 
surement in r2oo 
M200 

our best SIS fit to the shear data. The mass at r2oo is 17% lower 



7.05+1^2 X 10'"*^ 



With this model our total mass at r2oo is 
Q. The long -dashed curve represents 



in this case. The total X-Ray mass of Z3146 (see Sect. 3.5 1 is 
given by the short-dashed line. It is M200 = 9.57'014Mo. 
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Fig. 21. 7?-band image of Z3146. The image has a FoV of 
~ 16'40" X 16'40" (4.26 /i^g' Mpcx4.26 /i^,} Mpc), the Unes 
are the masks which were used to remove satellite tracks, fore- 
ground galaxies and bright halos from stars (see text for more 
details). North is up. East to the left. 



WEIGHT maps created by the data reduction pipeline (see 
lErben et al.|[2505j|Erben & Schirmer|[2003| for more details). 



sure that all objects are measured with the same photometric 
apertures. 

We used SExtractor in double image mode with the deep 
R band image as detection frame and the following param- 
eters: DETECT_THRESH=7, ANALYSIS_THRESH=7, and 
DETECT_MINAREA=3 (the higher detection threshold com- 
pared to the weak lensing analysis is a result of the seeing 
coiTection). All magnitudes are obtained using MAG_AUTO 
with PHOT_AUTOPARAMS= 1,3.5, as eUiptical apertures and 
a Kron radius of this size is best suited to our observations. 
In order to obtain clean catalogues with a minor fraction of de- 
fective detections like obvious stars/foreground galaxies, traces 
of asteroids and spurious detections in bright haloes of stars 
we masked such objects to remove them from the final cat- 



alogues. The image of Fig. 21 shows as example the original 
R-bwd image including all masked objects within a FoV of 
~ 16'40"x 16'40" (4.26 Mpc X 4.26/1^0 Mpc in our cos- 
mology). All masks were identical for the final B and V image, 
except for the individual satellite tracks. 
The total galaxy catalogue contains 2138 objects having a 
CLASS_STAR parameter < 0.95, MAG^UTO<99 in all bands 
(considering an E{B — V)- 0.126 mag, taken from the NEE|^ 
based on [Schlegel et aL] (1998 )) and a FLUX_RADIUS > 
3.2 pixels in B, V and R, respectively. In addition we used 

* http://nedwww.ipac.caltech.edu/ 



5.2.4. Analysis of the strong lensing features 

The results of the photometric and morphological investi- 
gations of all 4 arc candidates are summarised in Table|4] 
and Table|5] respectively. In this section we analyse the arc 
candidates using these informations. 

We can roughly estimate the strong lensing mass inside 
an Einstein ring at the position of the outer most arc Al 
(~ 27" ~ 120kpc). This mass can be estimated to be 
Md = 1.39!°i^^ X lO'^^Mo, where the main value is derived for 
(^source) - 0.79, the mean redshift value obtained in Sect. |5.1.2l 
The errors are calculated for Zsoui-ce = [2 x Zz3i46>2]. However, 
as we do not know either the redshift or the geometric align- 
ment of the source with respect to the lens, this procedure 
gives only a rough upper limit of the mass in the core. 
One of the most interesting questions is the possibility of 
finding multiple images of one single background source. 
Unfortunately we do not have spectra of the objects (see 
Sect. |5.2.r| which allow a secure identification of counter 
images. Hence we search for counter images in the following 
way: 

Counter images of arcs may not appear as elongated objects 
in the case of a folded arc system. In addition, they can differ 
in magnitudes due to the gravitational magnifying effect and 



can appear in unexpected locations ( Broadhurst et al. 2005a I, 
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which are not predictable without a precise model. Hence we 
have to restrict the identification of multiple lensed objects to 
investigations of the colour information (B - V), (B - R), and 
(V - R) only, as they are conserved by lensing. 
The search for multiple images was performed for all 4 
arc candidates independently in 2 steps: first, we searched 
the galaxy catalogue for objects with (a) coinciding colours 
(V - R), (B - R) and (B - V), and (b) lying in a radius of 3a.'0 
with respect to the cluster center position in the RBS. In a 
second step we discarded all objects being obvious cluster or 
foreground galaxies by visual inspection. 
arc candidate Al: In total we found two objects which might 
be counter images of candidate Al, denoted by CI, and C2, 



respectively (see Fig. 23 i. It is hard to judge whether those 
objects are counter images as they are hardly visible in the 
shallow WFPC2 image. However, C2 is located at a distance 
of ~ 30" with respect to the cluster center and hence we expect 
it to be much more sheared at this position if it had originated 
from the same object as Al. Additionally, the colours agree 
only within their large error bars (see Table|5]for the numbers). 
Hence we conclude that it is quite unlikely that this object is a 
counter image of Al. 

arc candidates A2 and A3: Both candidates show colour 
coincidences with each other and CI - C3. However, again 
the colour differences only agree within their error bars (see 
Table|5]for the numbers). The colour of A2 might be reddened 
to a certain amount by elliptical galaxies in its vicinity, 
nevertheless it is very unlikely that these objects have the same 
source. 

arc candidate A4: We did not find any counter image candi- 
dates for this object. 

Counter images also often occur nearby the central galaxy 
in the case of a not perfect alignment between the observer, 
the lens and the source. Hence they lie in the halo of a bright 
central object affecting their colour and/or are dramatically 
sheared up to a radial arc (see |Sand et aL| |2005| for some 
examples). Therefore we investigated the BCG using the HST 
image in more detail. 

A closer look at the WFPC2 exposure of the BCG reveals some 
knots in its very central part. To investigate these structures 
and to look for a possible radial arc we subtracted as a first 
step an elliptical model derived by fitting ellipses to isophotes 
of the BCG (done with the help of the IRAF tasks isophote 
and bmodel in the STSDAS package) as well as an artificial de 
Vaucouleur profile (task mkobjects in noao . artdata). 



Fig. 22 shows the central part of 73\A() with and without the 
subtracted elliptical isophote model of the BCG. The removal 
of the BCG reveals, apart from several clumps, an elongated 
substructure in the centre of the BCG along its major axis 
in the opposite direction to Al. However, we need deeper 
observations for identification of this object. At the cuiTent 
stage we can neither exclude the possibility of this object of 
being a radial arc or a filamentary structure common in cooling 
flow clusters. 



However, the fact that we did not find definitive counter 
images in our observations does not mean that there are none. 
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Fig. 23. R-bsnd image taken with WFI. Possible counter im- 
ages of the arc candidates are denoted with Cx. See Sect. |5.2.4| 
for the relations of the Cx to the arc candidates. 




Fig. 24. iV-R) vs. R colour-magnitude diagram. A distinct Red 
Sequence is visible, framed by two solid lines. The dashed lines 
show the limits used for identification of cluster members (see 
text for more details). 



Lensed sources can appear as very faint and thin arclets which 
are only visible in deep HST observations. Such arclets are 
therefore hard to find in ground based observations. Some 
prominent examples of lensing clusters such a large number 



of faint arc(lets) are e.g. A2218 (see |Soucail et al. 



2004 and 



references therein), A1689 dBroadhurst et al.[ |2005^i7A370 



(Bezecourt et al. 


1999) or CL0024-I-16 (|Broadhurst et al.| 


,2000„Kneib et al. 


|2003|). 
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Fig. 22. HST image of the innermost part of Z3146. In the right image we subtracted an elliptical model (based on fitting an 
ellipse to each isophote obtained by the IRAF tasks isophote and bmodel in the STSDAS package). Several distinct substructures 
are visible, but require deeper observations for their identification. 



arc candidate 


angular distance 


projected dist. 


length / 


l/w 




filter 




in Fig 


19 




to cc 


to cc [kpc] 


["] 


WFPC2 


B [mag] 


V [mag] 


R [mag] 


Al 




~ 27 


~ 117 


2'.' 3 


5.1 


24.49 ± 0.28 


23.87 + 0.06 


22.73 ± 0.01 


A2 




~ 20 


~ 87 


Z'O 


5.6 


23.89 ± 0.22 


23.07 + 0.04 


22.06 + 0.01 


A3 




~ 23 


~ 100 


v:a 


2.1 


24.20 ± 0.21 


23.68 ± 0.05 


22.71+0.01 


A4 




~ 26 


~ 113 


I'.'i 


1.9 


24.64 ± 0.22 


24.08 ± 0.05 


23.64 + 0.02 



Table 4. Table of the arc candidates shown in Fig 
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'cc' denotes the cluster centre, the length-to-width ratio l/w is calculated 
with the help of Ai, and A2 (the second order moments of the light distribution) and is measured on the WFPC2 images as it is 
not affected by atmospheric blurring. A value of = 0.75 was chosen for both, the DETECT_THRESH and ANALYSIS _THRESH 
due to the shallowness of the image. The magnitudes are instrumental WFI magnitudes in the Vega system. All objects are 
photometrised using SExtractor ' s MAG_AUTO. 



6. Investigations of the cluster light distribution 

In this section we present additional optical investigations on 
Z314 6 which are based on the WFI frames. 



6.1. Cluster member catalogue 

Independent of the previous analyses we have created different 
catalogues as we have different selection criteria for the further 
investigations. The lensing analysis focuses on background ob- 
jects, whereas the following investigations deal with the cluster 
members. We extracted a catalogue of cluster members in the 
following way: 



(a) From the galaxy catalogue created in Sect. 5.2.3 we made 
the colour-magnitude diagram (V - R) vs. R (see Fig. 24 1. In 



this plot we identify a Red Sequence (henceforth RS, marked 
by the two solid lines) which is used as the basis for the clus- 
ter member detection. The extraction of the Red Sequence was 



done by eye. As the RS galaxies belong to the ellipticals, which 
are the reddest ones in a galaxy cluster, we use the upper 
limit of the RS distribution as the natural colour border and 
assume all objects below the upper RS limit as cluster mem- 
bers. Additionally we skipped all objects with R > 23.5 mag as 
likely belonging to a background population, and objects with 
R < 18 mag as likely foreground systems. With these criteria 
we found in total 756 RS galaxies and 1478 cluster members. 

6.2. Galaxy distribution in Z3146 

To investigate the distribution of the RS members we created 
galaxy density maps in the following way: a blank image 
of about 4200 x 4200 pixels (corresponding to a FoV of 
~ 16'40" X 16'40", 4.26 h'^ Mpc x 4.26 h'^ Mpc) was 
created with pixel value "0" everywhere. At each position 
of the extracted Red Sequence galaxies (see Sect.|5]l the 
pixel value was changed to "1" and a subsequent Gaussian 
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object 


(V-R) 


A(V - R) 


(B-R) 


A(S - R) 


(B-V) 


A(S - V) 


Al 


1.14 


0.07 


1.76 


0.29 


0.62 


0.34 


A2 


1.01 


0.05 


1.83 


0.23 


0.82 


0.26 


A3 


0.97 


0.06 


1.50 


0.22 


0.52 


0.26 


A4 


0.44 


0.07 


1.00 


0.24 


0.55 


0.27 


CI 


1.15 


0.12 


2.37 


0.79 


1.22 


0.87 


C2 


1.02 


0.1 


2.26 


0.74 


1.24 


0.81 


C3 


0.89 


0.1 


2.74 


1.06 


1.85 


1.12 



Table 5. Photometiic properties of the arc candidates and their possible counter images. See Sect, 
discussions. 



5.2.4 



for more details and 




(y.'69 between the optical and the X-ray centre ([Schwope et al 



Fig. 25. Galaxy density plot of Z3146. The FOV (~ 16'40" x 
16' 40") is centred on the BCG, north is up, east to the left, 
I'.'O - 4.36 kpc. The image shows the distribution of 
the main Red Sequence galaxies shown in the upper plot of 
Fig 24 The contour levels coiTespond to 0.6,0.8, 1, ...,2.4 x 
lO^^alaxies/arcsec^. (The small dip in the contours between 
the cluster center and the bright star in the north is most likely 
an artifact from the masked and hence removed parts of the 
image, see Fig.[2T|) 



smoothing with cr - 241 pixels (corresponds to 250 /z^q kpc) 



leads to the image in Fig 25 



The galaxy density plot of the main Red Sequence (Fig 25 1 
shows one large peak centred on the main cluster with no 
distinct subclumps, except one small peak south of the cluster 
core. This is an indication that Z3146 is a relaxed cluster 
without any ongoing major merger event, which is confirmed 
by the massive cooUng flow found in previous investigations 
(Edgeetal. 1994 Fabian et al. 2002| l and our own results of 
~ 1600 Mq per year. In particular, the small distance of about 



|2000[ ) also confirms the calm character of this cluster. 



6.3. Light distribution / mass-to-iight ratio 

In order to obtain a mass-to-light ratio and creating a light 
distribution map we appUed the K-correction as a first step. 
We used the MatLab® script lum_f unc . m written by Eran 
Ofelj^ for this purpose. As iriput parameters we used the cor- 
responding WFI filter curveq^and template spectra provided 



Coleman 



by Stephen Gwyrj^ which are based on spectra by 
et al. ( jl980 ). As Red Sequence galaxies are mainly ellipticals 
we used E/SO spectra for them and Sbc templates for the re- 
maining. The resulting K-corrections (see Table[6]) were ap- 
plied to the galaxy catalogues of all galaxies in the FoV. To take 
the contamination resulting from non-cluster members into ac- 
count we created another catalogue of galaxies with the same 
criteria from a diff'erent region on the final WFI frames. This 
field is centered on a = 10*'22"59' and 6 = -h04°20'33.3", has 
the same size as the region used for creating the cluster mem- 
ber catalogue and has no distinct galaxy density peak. Hence 
we assume the galaxy population in this field to be dominated 
by field galaxies. Additionally this region is, in spite of the 
pointing off'set between the different observing programs (see 
Sect.|4|l, visible on all three observed bands. The galaxy counts 
in this field were binned in the same way as in the cluster field 
and subtracted from the corresponding bin of the cluster count. 

To calculate the total luminosity in the V band in solar units 
we assume the solar absolute magnitude to be My =4.82mag 
( |Cox[|20()0l l. We also assume the cluster members to follow the 
standard Schechter luminosity function ( Schechter 1976 1: 



(12) 



^{L)dL ^ f l^ — ^ e<i''>d{LIL*) 
In terms of magnitudes the Schechter function reads 
4>{M)dM = O.41n(lO)-0*-lO-°-4<^-^*>^''+'>-e[-"'""""*']^/M(13) 



7 


http : //wise-obs . tau. ac . il/~eran/matlab .html 




8 


http : //www. Is . eso . org/lasilla/sciops/2p2/ 


'E2i 


)2M/WFI/f liters/ 
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http : //orca . phys . uvic . ca/~gwyn/pz/specc/ 
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B 


V 


R 


units 


E/SO K-correction 
Sbc K-correction 
fitting range Mfj, 

f 
a 
M' 

x' 

Ltot 
F 


1.41 
0.84 

< -19 
1.814!|24x 10^ 

-1 136+°" 
-20.6!|i 
0.972 

1.12 X 10'^ 

< 1% 


0.85 
0.33 

< -19 

0.8!|^| X 10^ 

1 QO + I.30 

-21 74+3'" 

'^-3.01 

0.952 

0.98 X 10'"* 

< 1% 


0.32 
0.14 

< -19 
0.907^; X 10^ 

-1 24+"-" 
-22.53!];^^ 
0.973 

1.38 X lO'"* 

< 1% 


[mag] 
[mag] 
[mag] 

[# galaxies deg"-] 
[mag] 



Table 6. Details of the catalogues, K-coiTections values and the Schechter function parameters 0*, a and M* for the three WFI 
filters. F is the fraction of light which is missing due to the limiting magnitudes M\i^ of the catalogue given in Table[3] For the 
R band we used the limiting magnitude for region A (see Table[3]). See text for more details. The Luminosity distance is 1.5Gpc 
(distance modulus m - M = 40.88 at z = 0.2906). 




The total luminosity L,,,, (see Table[6]l can now be obtained by 
integrating Eq.[T2| 



L,o, - 0*LT(2 - a) 



Fig. 26. V band light distribution map. The FoV is ~ 16'40" x 
16'40", North is up. East to the Left. The contours correspond 
to values 2 x 10^ 2.5 x 10^ 3 x 10^ 7.5 x 10^ 8 x 10^ Lq. The 
light distribution shows one single peak without any distinct 
substructures. The position difference between the cluster cen- 
ter and the peak of the distribution is about ~ 29", which might 
be a result from the masked parts of the image, see Fig.[2T|i 



To obtain the parameters L*, (jf and a we applied a fit to 
Eqjisjusing the MatLab® Fitting Toolbox at a 95% confidence 
level. The best fitting parameters including the ;\f^-value for 
the goodness of the fit and the fitting range for the luminosity 
[inf, L/„] are given in Table|6] where L/„ is the completeness 
luminosity. 



(14) 



The integration of the Schechter function down to a luminosity 
Liim is equal to 



Liim - ^*L*r{2 - a,Lii,„/L*). 



(15) 



As the Schechter function is not applicable to very faint lumi- 
nosities we calculate the fraction of light we miss due to obser- 
vational effects 



1 



r(2-a,Lf,/L*) 



(16) 



r(2-a,L,,„,/L*) 

Lfit being the luminosity (corresponding to M/„ in Table|6|l 
down to which our catalogues are complete and L;„„ is the lim- 
iting luminosity corresponding to the limiting magnitude given 
in Table|3] 

Fig|26] shows the corresponding light distribution of the Red 
Sequence galaxies in V. Again, only one single, very distinct 
peak centred on the BCG is visible. In addition the distribution 
is very smooth and shows no substructures, confirming the re- 
laxed state of Z3146. 



In Sect. 5.1.2 we fitted the mass obtained by the weak lens- 



ing method to an NEW profile with the following best-fit pa- 
rameters obtaining r2oo - 1661^32" /i^dkpc and c = 3.6;|;2 4 
as best fit values. This leads to M200 = 7.04+^']^ x IO'^^Mq. 
Using these values we find a mass-to-light ratio within r2oo in 
the V band of M/Ly ~ 156+^^'' Mq/Lq. This value is in agree- 
ment with e.g. Hradecky et al. ( |2000[ ), who give a median value 
of M/Lv ~ 138Mo/Lq for eight clusters within a radius of 
~ I.38/17IMPC. 



7. Discussion and summary 

We presented a combined investigation of optical and X-ray ob- 
servations of the prominent galaxy cluster Z3146. This cluster 
seems to be in a relaxed state, which is confirmed by 
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- the absence of large substructures in the galaxy density plot 



and the light distribution map (see Fig. 25 and Fig. 26 1, 

- the regular shape of the cluster in the X-ray image and the 
temperature map (see Fig.[T]and[3] respectively) 

- the massive nominal cooling flow of ~ 1600Moyr"' 

- the good coincidence of the optical, the X-ray, and the weak 
lensing centre (each of the order of a few arcseconds), and 

- the regular shape of the weak lensing mass reconstruction. 

Further optical investigations on the cluster also revealed 
four gravitational arc candidates and a mass-to-light ratio of 
MILv ~ 156!if h-l^Me/L^ at rzoo = l^bVj^ /z.-jkpc. 

We also determined the mass of this cluster with two inde- 
pendent methods, weak lensing and X-ray measurements. Both 
data sets. X-ray and lensing, were used to establish best fits to 



the commonly used NFW model. Fig. 27 shows the compari- 
son of the confidence levels for these data: 
Fig.^7\a) - Xray-data: Confidence level from a NFW fit to 
the mass profile derived from X-ray data according to Eq. |4] 
The contour levels are 2.31, 6.25, 1 1.90 corresponding to con- 
fidence levels of 68.3% (Icr), 95.4% (2cr), 99.73% (3cr). 
Fig.^7\b) - Lensing data: The contours are at 2AX =2.30, 6.17, 
9.21 coiTesponding to confidence levels of 63.8%, 90%, 95.4% 
and 99% if we assume Gaussian statistics. We varied the galaxy 
sample of the lensing fit to investigate the dependence of the 
result on this parameter. On the one hand, lowering the maxi- 
mum radius to which galaxies enter the calculations to 8'.0 from 
the galaxy centre (we have full data coverage around the clus- 
ter up to this radius; see Fig.jsjl leads to r2oo = 1610 /i^^^kpc; 
c = 4.01. On the other hand, not considering the inner parts 
of the cluster and using only galaxies with a distance larger 
than 2'.0 we obtain r2oo = 1748 /i^^kpc; c = 1.0. We note that 
r2oo is reasonably well constrained and that the concentration 
c mainly depends on the details near the cluster core. This be- 
haviour coiTesponds to the shape of our contours and is typical 
for NFW profile fits in weak lensing studies; see e. g. |Clowe| 
& Schneider] p002| and [DieSchetirLl ( f2005| ). The parameter 
ranges in c and r2oo imply an uncertainty of the total cluster 
mass of 10% - 20% (considering radii of 1 - 2 hj^Mpc; see 



also Fig. 18 1. 

Fig.^7\c) - Direct Comparison: The 'x' marks our best fit lens- 
ing value of r2oo = 1661 /i^q kpc (r^ - 460.1 kpc) and 
c - 3.61, which lies in the vicinity of the 3cr X-ray model. 
The triangle corresponds the best NFW fit to the X-Ray d ata: 
3.58 (r, 



'"200 - t»4V '170 i^pc, c - (rs - 016.6 /i^q kpc, see Sect. 3.5 
for more details). This best fit value is located within the Icr 
contour of the lensing model. Hence both models are in excel- 
lent agreement. 

A direct comparison of the mass profiles and the ratio be- 
tween Mx-rav/Mwi is given in Fig. 28 which shows that the best 
fit models agree within ~ 20%. 

Comparing the strong and weak lensing masses it seems that 
they disagree. At the radius of the outermost arc at r ~ 27" 
(~ 120 kpc) we obtain a mass within an Einstein ring of 
about = l-39+[;;^^ x lO'^ h-^M^ 



for the strong lensing 
measurement. The NFW profile of the weak lensing fit gives 
Mwi - 2.351? 5? ^ ^^'"^ ^70 position, which is. 






3ICC- 



Fig. 27. Comparison of confidence levels from NFW fits: Panel 
(a) shows the fit to the X-ray data. Panel (b) the fit to the shear 
data. A direct comparison with the corresponding best fit val- 
ues is presented in Panel (c), where the lensing model is repre- 
sented by the solid lines (best fit marked by 'x'), and the X-ray 
model by the dash-dotted lines (best fit values correspond to 
the triangle). See Sect.|7]for more details. 

the best case assuming, roughly half of Mji only. However, due 
to the large uncertainties in the strong lensing mass determi- 
nation (unknown redshifts, unknown lensing geometry...) we 
assume this mass only to be a rough upper value. Hence this 
discrepancy is likely an artifact of the large numbers of uncer- 
tainties in the determination of the strong lensing mass. 

Especially in relaxed clusters, the mass estimates obtained 
from weak lensing and X-ray mass methods usually seem 

|Wu et al.| [T998|. Recent 



to agree very well (Allen 1998 



observations of cooling flow clusters derived from Chandra 



and/or XMM-Newton confirm these results (see e.g. Allen 
[eFaH [2OO2I [Cypriano et"aL| [2005] ). We find a temperature of 
5.9 + 0.1 keV in Z3146, in agreement with the assumption of 



Cypriano et al. ( 2004) that clusters having an ICM temperature 



< 8.0 keV are in a relaxed state. In particular relaxed clusters 
are interesting for cosmological studies as their mass content 
tends to take a spherically symmetric shape, which is the usual 
assumption in theoretical approaches. Hence a large sample 
of such systems is a useful probe to verify whether the mass 
density of galaxy clusters follows an NFW profile ( [Navarro 
et al. 1996 1997 1, or whether a different profile like the 
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Both show good coincidences within Icr errors (see Fig. 27 1. In 



the lower panel the ratio of Mx-ra\IM„i is shown. Note, that the 
fit of the X-ray data to a NFW profile starts at < 230 /!^(Jkpa 
The dash-dotted line represents r2oo 
for more details. 



1661/!7o'kpc. See Sect.]?] 



Burkert ( [Burkert| |2000| i, the Moore ( |Moore et aL]|1999| l or the 



non-extensive profile ( [Leubner 2005 [ Kronberger et al.j 2006| l 
is a suitable description. 
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